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Abstract 



We consider an integrable generalization of the sine-Gordon (sG) equation that 
was earlier derived by one of the authors using bi-Hamiltonian methods. This 
l_H equation is related to the sG equation in the same way that the Camassa-Holm 

^0 equation is related to the KdV equation. In this paper we: (a) Derive a Lax pair. 

£H (b) Use the Lax pair to solve the initial value problem on the line, (c) Analyze 

solitons. (d) Show that the generalized sG and sG equations are related by a 
Liouville transformation, (e) Derive conservation laws, (f) Analyze traveling- 
wave solutions. 
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1 Introduction 

ON 

We consider the following integrable generalization of the sine-Gordon equation, 
which was derived in [6] using bi-Hamiltonian methods: 

.it (1.1) utx = (1 + i/fl£)sin(u), t > 0, 

where u(x, t) is a scalar-valued function and v S R is a parameter. Equation 



1.1) 



is related to the sine-Gordon equation (i.e. to the equation obtained from ( |1.1 ) by 
letting v = 0) in the same way that the Camassa-Hom (CH) equation (see OUT]) is 
related to the Korteweg-de Vries (KdV) equation. Actually, there exist even deeper 



analogies between (1.1) and the CH equation. Indeed, recall that for a particular 
class of initial conditions it is possible to use a Liouville transformation to map CH 
to an equation in the KdV hierarchy. Similarly, it will be shown here that for v < 0, 
equation ( |1.1[ ) is related to the sine-Gordon equation via an appropriate Liouville type 



transformation (see Proposition 5.1). Furthermore, one of the distinguished features 



of the CH equation is that it possesses certain non-smooth traveling-wave solutions 



called peakons. Similarly, it will be shown here that equation (1.1), in addition to 
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smooth traveling- wave solutions (see Proposition |7.1[ ), also possesses certain non- 
smooth traveling- wave solutions called cuspons (see Proposition 7.2). 

Taking into consideration the extensive interest in the literature for the CH equa- 
tion, as well as the above similarities between ( |1.1[ ) and CH, it is natural to expect 
that equation (1.1) may also be of some interest. 

This paper is organized as follows: Section [2] presents an algorithmic derivation of 
a Lax pair for equation (1.1), starting from the knowledge of the recursion operator 
(this approach might be useful also for other integrable equations generated via bi- 
Hamiltonian methods, since for such equations the associated recursion operator is 
known) . Section [3] presents the solution of the initial- value problem of equation (1.1) 
with decaying initial conditions in the case that v < 0. The more interesting case of 
v > remains open. Solitons and conservation laws are discussed in sections [4] and 
[6] The explicit Liouville transformation relating (1.1) with the sine-Gordon equation 
is presented in section [5] Smooth, as well as non-smooth traveling- wave solutions are 
derived in section 

Taking into consideration that ( 1.1 ) is related to the sG equation with an explicit 
transformation (for v < 0), it follows that section [3] also presents the solution of 
the initial-value problem of the sG equation. The solution of this problem was first 
presented in the classical work pQ, however the method presented in section [3] is based 
on the 'modern' Riemann-Hilbert approach and is also influenced by the method of [7] 
(so that it is not necessary to formulate separately the time evolution of the spectral 
functions, see [TO]). 

After this work was completed, the references j2[ [16] (where this equation is stud- 
ied in connection with pseudospherical surfaces) and [IT] (where the transformation 
between ( 1.1 ) and the sine-Gordon equation is presented) were brought to the authors 
attention. 



2 A Lax pair 
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Equation (1.1) is the compatibility condition of the following Lax pair: 

J ip x + iAo 3 V> = W%ip, 



(2.1) 



where if)(x,t,\) is a 2 x 2-matrix valued function, AeC = CUooisa spectral 
parameter, and Wj(x, t, A), j = 1, 2, are defined by 



(2.2) 



Wi =u x ( -01 + X^a 2 ) , 



W 2 =— [sin(n)o 2 + (cos(n) - l)o 3 ] + 
4A 



IV 



cos urn 



x 



sin(it) 



01 



+ A 



V 3 / 2 COs(u)u x (T2 — ^(COS('U) — l)o 3 
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2.1 Derivation of Lax pair 



Before deriving the Lax pair (2.1 ), we first recall the case of the sine-Gordon equation 
(2.3) q t = sin(d- 1 q). 



Equation (2.3) admits the isospectral problem (x-part of a Lax pair) 

i 

(2.4) tp x = Lip, L = -i\a 3 + -qa%, 

where ip(x, X) is a vector valued eigenfunction and A € C is a spectral parameter, as 
well as the recursion operator R defined by 

R = e 2 d-\ e 2 = di + d x qd- l q d x . 



Differentiating equation (2.4) with respect to q in the direction w (where ip and A are 



considered functionals of q), we find 



d x Tp'(q)[w] = Lt/}'(q)[w] - i / G x wdxa 3 i/j + -wa^, 



(2.5) 

where G\ denotes the gradient of A defined by 

A'(g)[u;] = J G\wdx for all w. 
Let tp* be a solution of the adjoint equation 

(2.6) r x = -LV, 

where L A = L denotes the adjoint of L with respect to the inner product 
(f,9) = j fgdx, i.e. (f,Lg) = (L A f,g). 



Multiplying (|2.5l) by ip* T from the left and integrating the resulting equation with 



respect to x, we find 

(2.7) / ip* T d x ^'{u)[w]dx = I i;* T L^'{u)[w}dx 



+ 



G\[-i I tp* T a 3 ipdx ) + -i/j* 1 aiip 



i.*T. 



wdx. 



The first two terms of (2.7) cancel in view of (2.6). Thus, since w is arbitrary, 



(2.8) 



G x 



2 J ip* T asipdx 



From the bi-Hamiltonian theory, we expect the gradient G\ to satisfy an eigenvalue 
equation of the form [Hj 



(2.9) 



R A G X = K\)G> 



where \i = /i(A) is an eigenvalue and R = d x 1 02 denotes the adjoint of R. In fact, 
letting ip = ip* T a\ip, we find that ip satisfies the equation 



/ n(X)tp - ip x 
V Q 



M(A) 



-4A 2 , 



which simplifies to (2.9) (up to the irrelevant constant in the denominator of ( |2,8[ )), 
Equation (2.9) shows how the eigenfunctions of the x-part are related to the 
recursion operator for the sG equation. The eigenfunction <p = ip* T a\ip is sometimes 
referred to as a 'squared eigenfunction' since it is bilinear in ip* and ip. 



We now proceed to the derivation of the Lax pair (2.1). Equation (1.1) possesses 
the recursion operator [6] 



R — 020\ 



where 



dl + d x u x d x 1 u x d x . 



In order to find a Lax pair of the generalized sine-Gordon equation (1.1) we have 
to trace the steps used for the derivation of (2.9) backwards. Namely, knowing the 
recursion operator R = 020 \ X , the goal is to determine the corresponding x-part. 
Thus, we seek a n x- part such that the gradient G\ of the spectral parameter A of this 
x-part satisfies (2.9) with R A = 0^ 1 02- We let q = u x and make the ansatz 

ip x = Lip, L = -i\a 3 + f(q, A)oi + g(q, A)cr 2 , 

where / and g are scalar-valued functionals of q and of the x-derivatives of q which 
depend on A. The analog of equation (2.5) is 

(2.10) d x ip'(q)[w] =Lip'(q)[w] -i [ G x wdxa 3 ip + 



G\wdx&\ip 



+ 



dg_ 
dX 



G\wdxa2ijj + (Dfiv)aiip + (D g w)a2ip, 



where the operators Df and D g are defined by 



f dq dq x x dq xx 



8i 81 + 



D _ dg_ dg_ g _d§ 
9 dq dq x x dq x 

*T 



-dl + 



Letting ip* be a solution of (2.6), multiplying (2.10) by ip* from the left and inte- 
grating the resulting equation with respect to x, we find the following equation: 



ip* T d x ip' (u) [w]dx 
(2.11) + 



ip* Lip' (u)[w]dx 



G\[-i I ip* T a 3 ipdx + / ip 



,.*rdf 
dX 



o\ipdx 



,*T d 9 , , 
W -7y K °2Wdx 



+ DfiiP* 1 a x iP) + D^iP* 1 a 2 iP) 



wdx. 
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The first two terms of this equation cancel and hence G\ is given by 



T dR. 



Letting 



f(q,X) = a(\)q, g(q, A) = b(X)q, 

where a(A) and 6(A) are complex-valued functions of A, we find that the function (p 
defined by 

satisfies the equation 
(2.12) 
whenever 

(2.13) M (A) 



( + V-Pxx) ~ <fxx 

q 



4A 2 



1-4^A 2 ^ 



a = % -\/\ + 46 2 -4A 2 i/. 



Rearranging (2.12), we deduce that G\ satisfies (2.9) whenever fi,a,b satisfy (2.13). 
Choosing for simplicity a = i/2 and b = we find 



i>x = Lip, 



L = -i\a 3 + -qai + \y/vqo2- 



Having obtained the x-part of the Lax pair (2.1), in order to find the corresponding 
t-part, we make the ansatz 

i 

i> t = Mi>, M = A J '(A iC Ti + Bja 2 + Cja 3 ), 
j'=-i 

where Aj,Bj, Cj are scalar- valued functionals of q independent of A. Identifying terms 
of 0(X J ), j = —1, 0, 1, 2, in the compatibility equation L t — M X + [L, M] = and using 
that qt = (l + i'Q 2 .) sin(9~ 1 g), long but straighforward computations lead to the i-part 
in (1211). 



3 Spectral analysis on the line 

The Riemann-Hilbert (RH) formalism for integrating a nonlinear evolution equation 
is based on the construction of eigenfunctions of the associated Lax pair which can be 
joined together to a bounded and sectionally analytic function on the Riemann sphere 
of the spectral parameter A 6 C. In order to define eigenfunctions which possess the 
asymptotics 

/j = I + 0(l/\), X-^oo, 



we will use a different representation of the Lax pair (2.1). This representation 



involves an eigenfunction (f>, which is related to ip via the gauge transformation ijj 
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g<j). The form of g is such that it (a) diagonalizes the highest order terms in A as 
A — > oo and (b) preserves the natural symmetry properties of the Lax pair. 

The spectral analysis of the Lax pair satisfied by <j> is standard: We introduce two 
particular eigenfunctions via integration from x = ±00, and use these eigenfunctions 
to formulate a 2 x 2-matrix RH problem with jump across the real axis. The jump 
matrix is given in terms of two spectral functions a(A) and b(X), which are defined 
in terms of the initial data u$(x) = u{x, 0). The possible poles of a(A) give rise to 
singularities of the RH problem which correspond to solitons. 

The solution u(x,t) can be recovered from the solution of the RH problem, so 
that this approach provides the solution of the Cauchy problem for equation 
However, the RH problem is naturally formulated in terms of variables (y, t) rather 
than (x, t) where y is a new variable. Therefore, we will only obtain the solution 
u(x, t) in parametric form. This type of parametric representation arises also in the 
analysis of the Camassa-Holm and Degasperis-Procesi equations (cf. [U [12] in the 
case of the Camassa-Holm equation). 

In the rest of the paper we assume that v < and then for simplicity let v = — 1. 



3.1 Symmetries 



Assuming that v = —1, the Lax pair (2.1) can be written as 
(3.1) 
where 



ipx = Lip, 
ipt = Mil), 



L -i\cr:\ + \\\, M ! A - 4A ) • II. 



and W\, W2 are given by (2.2) with y/v replaced by i. Since L and M are trace-less 
and satisfy the symmetry relations 

£t(A) = -L{\), M\\) = — M(A), 

where denotes the complex-conjugate transpose of a matrix A, it follows that the 
functions 

det(ip(x,t, X)) and ip'(x,t, X)ip(x, t, X) 

are independent of x and t. We deduce that ip satisfies 

(3.2) det(i/}(x, t, A)) = 1, $(x, t, X) = ^ _1 (x, t, A), 

for all (x,t) provided that these equalities hold at only one point (xo,to). The eigen- 
functions introduced in the sequel will indeed have this property. 



3.2 Eigenfunctions 



Suppose that u(x,t) is a sufficiently smooth real-valued solution of (1.1) and suppose 
that cos(u(x,t)) — 1 has decay as x — > ±00 for all t > 0. Then the form of equation 
(1.1) implies that J R sin(u(x,t))dx = for each t. Define m(x,t) by 

(3.3) m(x, t) = 1 + u x {x, t). 
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The gauge transformation 

(3.4) rP(x,t,\) = g(x,t)<f>(x,t,X), 

where 



(3.5) 



g(x,t) 



1 + y/m 
2,/rh 



1 i 

■ \—\fm 



1-ym 
1 



transforms the Lax pair (2.1) into 
(3.6) 



x + i\p x (T3(f> = Vl(j), 

1 \-*- V 



y </> t + i(\pt- 4xjcr3 

where Vi(a;,f) and V2(x,t, A) are defined by 



T _ 1/ / ^XX 

vi =— Ux H I cii 



i cos(u) — ti x sin(u) u x cos(u) + sin(it) 

7Y " 3 + * T7=^ ^3 + 1 TT=~\ a2 

4A AJmX AJmX 



V 2 ■■ 

% f ^'XX I / \ 

- n\ U x-\ COS(li)CTi 

and p(x, i) is a real-valued function defined by 

(3.7) p(x,t) = x -t+ [ (y/m(x',t) - l)dx. 

J — oo 

The conservation law 
(3.8) 



(v / m) i + (cos(u)y/m) = 



implies that pt = — cos(u)y / m. Note that definition (3.3) of m implies that g is free 
of singularities also at points where u x = 0. 



The form (3.5) of g is motivated by the fact that it diagonalizes the terms of O(A) 
of the Lax pair (2.1) and that it satisfies 



(3.9) 



det(g(x,t)) = 1, g ] {x,t)=g 1 (x,t). 



The relations (3.9) ensure that the gauge transformation (p.4| preserves the two 



properties in (3.2), i.e. 

(3.10) dek(<j>(x,t,\)) = l, t, A) = <p~ x {x, t, A). 

We have defined V\ and V2 so that 

Vi(x,t) ->• 0, V 2 (M, A) -)• 0, a; -)• ±00. 
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The form of the Lax pair (3.6) is convenient for the definition of eigenfunctions 



which are well-behaved near A = oo. Introducing an eigenfunction by 



we find that the Lax pair (3.6) can be written as 
(3.11) 



We define two eigenfunctions n± of (3.11 ) as the solutions of the following two Volterra 
integral equations 



ti+(x,t,\) = I- / e iXMx '' t) - p( - x ' t))a W 1 (x',t) f i + (x',t,X)dx', 

J X 

H-(x,t,\) = I+ [ X e lX(p(x '' t) - p{x ' t))a W 1 {x',t)^(x',t,X)dx', 



(3.12) 



where (T3 acts on a 2 x 2 matrix A by &3A = [03, A]. The second columns of these 
equations involve the exponential e 2lX ^ x ' >*)~p( x >*)) . Since p(x, t) is an increasing func- 
tion of x for any fixed t, we deduce that the second columns of fi + and /U_ are bounded 
and analytic for A in the upper and lower halves of the complex A-plane, respectively. 
We will use the superscripts (+) and (— ) to denote these boundedness properties. 
Similar considerations are valid for the first columns and hence 



The Lax pair (3.11) is of 'standard form' as A — > 00, i.e. the highest-order terms 



of 0(A) are diagonal and appear in the exponents of (3.12), whereas the next order 
terms of 0(1) are off-diagonal. Substitution into (3.11) of the expansion 



/x = D + ^ + ~ H , A^oo, 

A A z 

shows that D is a diagonal matrix independent of x and t. In particular, n + and /x_ 
have the following asymptotics near A = 00: 

A 4 - ? M+ 



J + 0(1/A), A 00, ImA>0, 
J + 0(1/A), A ^00, ImA<0. 
Note that //+ and fi- are nonsingular near A = 0, because the integral equations 



f 4 + ? P- 



(3.12) involve only the function V\ and not V2. 



3.3 Spectral functions 

We define the spectral function s(A) by 

(3.13) A*+CM,A) = / u-(x,t,A)e- i(Ap(3; '' ) -A ) " 3 s(A), 



ImA = 0. 
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Evaluation of (3.13) at t = and x — > — oo gives 



(3.14) s(X) = I 



3 iXp ( x >V &3 Vi(x,0)n+(x,0,X)dx, ImA = 0. 



We deduce from (3.10) that fi± have the properties 

(3-15) 

det(/i . (x, I, A)) = 1, n±{x,t,\)n = M±(x,i,A) 2 2, H±{x,t,X) 2 i = -//±(x,i, A)i 2 . 

This implies that 

dets(A) = 1 

and that there exist functions a(A) and 6(A) such that 



S (A) 



a(A) 6(A) 
-6(A) a(A) 



From the explicit expression (3.14) for s(A) and the fact that the second column of 



fi + is defined and analytic in ImA > 0, we find that a(A) has an analytic continuation 
to the upper half-plane. 

3.4 Residue conditions 

We assume that a(A) has N simple zeros {\j}j = i in the upper half-plane. These 
eigenvalues are either purely imaginary or arise as complex conjugate pairs (A, —A). 



The second column of equation (3.13) is 

,(+)_„„(-),,,„.(+) „-2i(Ap-A) 



(3.16) 



ImA = 0. 



Applying det ( \i 



,(+) 



to this equation and recalling (3.15), we find 



det A) ^\x,t,X)) =a(X), 



Im A > 0, 



where we have used that both sides are well-defined and analytic in the upper half- 
plane to extend the relation to ImA > 0. Hence, if a(Xj) = 0, then fi_ (x, t, Xj) and 

fi+ (x,t, Xj) are linearly dependent vectors for each x and t. It follows that there 
exist constants bj such that 

(3.17) n^\x,t,Xj) = b j e 2i(X3P ~^\ ( +\x,t,X j ), x G R, t > 0. 



Recalling the symmetries in (3.15), the complex conjugate of (3.17) is 



fi_ \x, t, Xj) = —bjt 



-2i(X jP -^) (_) 



j /i\ (x, t, Xj), iet,t>0. 
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Consequently, the residues of fi_ /a and fi_ /a at Xj and Xj are 

Res^%^ = C^^-^^^X,), 
^ a(X) 

He / L-W,A) = _^ e -»(^)-^) M H (a;)tiXi)) 
a j a (A) 

where d = 4r and C 



dA J "(A.,) 



3.5 Riemann-Hilbert problem 



We now describe how the solution to the initial- value problem for equation (1.1) on 
the line can be expressed in terms of the solution of a 2 x 2-matrix Riemann-Hilbert 
problem. Relation ( 3.13| ) can be rewritten in the form of the RH problem 

M-{x,t, A) = M + (x,t,X)J(x,t,X), ImA = 0, 

where the matrices M_, M+, J are defined by 

(3.18) M + = ^^y,/4 +) ^ , ImA>0; 

M_ = t— \ , ImA < 0; 



1 



a(A) 



(3.19) J = I tjts t , y > , ImA = 0. 

a ( A ) a(A)a(A) 

The contour for this RH problem is the real axis. 

However, since the jump matrix J, in addition to the spectral functions a(A) and 
6(A), also depends on the function p(x, t), this RH problem cannot be formulated in 
terms of the initial data alone. We overcome this problem by changing variables 

(x,t) h-> (y,t), y=p(x,t). 

Since J only depends on x through y = p(x, t) we can formulate a family of RH 
problem parametrized by (y, t). From the solution M(y, t, A) of this RH problem, we 
recover u(x, t) in parametric form. 
Defining 

/ i _M e - 2i (^-^)\ 

(3.20) jC>(y,t,A)= Rawa,-A) ^ i ' ImA = °< 



,bW e 2i(\y-^) 

a ( X ) a(X)a(X) 



we have the following result. 
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Theorem 3.1 Let uq(x), be such that cos(uo(x)) — 1 has sufficient smoothness 

and decay as x — > ±00 f/or details see (3.35) below) and such that f R s'm(uo)dx = 0. 
Define Vio(x), mo(x), and po(x) by 



m 



V w = ~ ( uox + V ^ J 0"i, "T,o = l + iiox' p (x) = x + (^m (x r ) - l)dx'. 



Let n+(x,0, A) and fJ>-(x, 0, A) 6e i/je unique solutions of the Volterra linear integral 
equations (3.12) evaluated at t = 0, mf/i Vi(x,0) = Vio(x) and p(x,0) = po(x). 
Define {a{X),b(X),C j } by 



s(\)] 2 , S (X) = L- e iX ^ a W 10 (x)^ + (x,0,X)dx, 



1,...,N, 



and 

(3.22) [M^A,)]! = aCAjOCje 2 *^^^^^^, A,)] 2 , 3 = 

where [A]i and [.A] 2 denote the first and second columns of a 2 x 2-matrix A. We 
assume that a(A) /ias JV simple zeros {Xj}jLi in the upper half-plane. 
Then 

• a(A) is defined for ImX > 0, analytic in ImX > ; and a(A) = 1 + 0(1/A) ; 
A — > 00. 

• 6(A) is defined for A £ E and 6(A) = 0(1/A), A — > 00. 

• |a(A)| 2 + |6(A)| 2 = 1, AGE. 

Suppose there exists a unique solution u(x,t) of equation (1.1) with initial data 
uq{x) such that cos(u(x,t)) — 1 has sufficient smoothness and decay as x —> ±00 for 
each t > 0. Then u(x,t) is given in parametric form by 

(3.23a) u(X(y,t),t) = 21m (j* a(y' ,t)djA + 2nvr, 

where 2mr = lim a; _ s ._ 00 uq(x), the function X(y,t) is defined by 

(3.23b) X(y,t) = y + t + [ ( y/l - 4Im(a(y', t)) 2 - l) dy' , 

J— 00 



i/te function a(y,t) is the unique solution of the Ricatti equation 

(3.23c) a y = a 2 — 4i ( lim XMi2(y,t, X)) a — - , lim a(y, t) 

\ .\ • \ y 4 . j- . 



y— »±oo 



1 

2' 



which satisfies 



Hy,t)\ = -, yeR, t>0, 



and M(y, t, A) zs the unique solution of the following RH problem: 
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M(y,t,X) 



M+(y,t,X), 
M-(y,t,X), 



ImA > 0, 
ImA < 0, 



is sectionally meromorphic. 



M-(y,t,X) = M + (y,t,X)j( y \y,t,X) for X G R, where jM is defined in 
M(y, t, A) has the asymptotic behavior 

1 



(3.24) 



M(y,t,X) = I + 



X — > oo. 



The first column of M + has simple poles at X = Xj, j = 1, . . . , N, and the 
second column of M- has simple poles at X = Xj, j = 1, . . . ,N. The associated 
residues are given by 

(3.25a) Res[M(y,t,A)]i = C^^^Mfo^Jh, 

(3.25b) Res[M(y,t,X)} 2 = -Cje' 2 ^'^ [M(y, t, A,-)]i. 



Proof. It remains to prove (3.23). By substituting the expansion 



fi(x,t, A) = / + \ ' + - A V + ° ( A3 ) ' A - x ■ 



into the x-part of (3.11), we find by considering the terms of 0(1) that 



(3.26) 



u x (x,t) + 



m(x, t) 



Define Q(y,t) by Q y (p(x,t),t) = —A/i^ix^t). Then, by construction of the RH 
problem, 



(3.27) 



Q y (y,t) = -4 lim (XM l2 {y,t,X)) 



X— ¥oo 



Equation (3.26 ) is conveniently expressed in terms of y = p(x, t) and Q. Indeed, using 
that d x = \fmdy and m = 1/(1 — Uy), we find 



(3.28) 



Qy — I Uy -\- 



U 



till 



Integration of this equation yields 
(3.29) Q = —u — arcsin^y), 

where we have fixed the integration constant by assuming that 

Q(y,t) + u(y,t) 0, y -> ±oo. 



12 



Note that 



-1 < u 



u , 



y 



< 1, 



so that Q is well-defined by (|3.29|). We rewrite (|3.29|) as 
(3.30) 

Defining a(y,t) by 



u„ 



sin(Q + u). 



(3.31) 



a 



AQ+u) 



equation (3.30) becomes the Ricatti equation 
(3.32) 



a y = a 2 + iQ y a - -, a(±oo,t) = — -. 



By ( |3.30[ ) and ( |3.3l| ), 

(3.33) Uy = 2Im(a). 

On the other hand, y — > x — t as x — > — oo, and so 

(3.34) 



x = y + t + 



1 - vl - 1) dy'. 



The parametric formula (3.23) now follows from equations (3.27), (3.32), (3.33), and 
(pT34|). □ 



The jump matrix defined in (3.20) satisfies the symmetry relation jj^ (A) = 
Jgj' (A). This implies that there exists a 'vanishing lemma' for the RH problem. This 
guarantees existence in Hi(M) [9]. Thus the initial conditions must be chosen in an 
appropriate function space such that 



(3.35) 



p (x) - x £ Loo, Vio e L 1: a,66Fi. 



4 Soliton solutions 



The solitons correspond to spectral data |a( A), MX ), Cj} for which 6(A) vanishes iden- 
tically. In this case the jump matrix j( y > in (3.20) is the identity matrix and the RH 



problem of Theorem 3.1 consists of finding a meromorphic function M{y,t,X) satis- 



fying (3.24) and the residue conditions (3.25). In what follows we will suppress the 



(y, independence of M and write M(A) for M(y, t, A). 



Let 9j = Xjy 



t 

4A, 



From (3.24) and (3.25a) we find 



(4.1) 



[M(A)] 



N 



+ E 



C s ^[M(Xj)]2 
A - A,- 
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It follows from (3.15) that M(A) respects the symmetries 
(4.2) 



Mn (A) = M 22 (A), M 21 (A) = -M i2 (A). 



Moreover, if (4.2) holds, then only one of the two residue conditions ( 3.25a )-( 3.25b) 



needs to be verified since the other condition is a consequence of symmetry. In view 



of (4.2), equation (4.1) can be written as 



(4.3) 



M 22 (A) 
-M 12 (A) ( 



C,e 2 ^ (M 12 (Xj) 
< x A - Xj \M 22 (Xj) 



Evaluation at A n yields 
(4.4) 



M 22 (A r 



-M 12 (A n ; 



A | y C 3 e™> (M 12 {\ 3 ) 
0/ ~^ A n — Xj \M 22 (Xj) 



N. 



Solving this algebraic system for Mi 2 (Aj) and M 22 (Aj), j = 1, . . . , N, and substituting 
the resulting expression into (4.1) and (14. 31) gives explicit expressions for the entries 



of M(A). The iV-soliton solution u(x,t) is then obtained from (3.23). 



4.1 One-solitons 

The one-solitons arise when a(A) has one purely imaginary zero Ai. Letting N = 1 



in (4.4) and taking the complex conjugate of the second row, we find 



Ci e 2idl 

M 22 (Ai) = 1 + -=— — — Mi 2 (Ai), 
Ai — Ai 



-Mi 2 (Ai 



Ai - Ai 



M 22 (Ai). 



Solving these equations for Mi 2 (Ai) and M 22 (Ai), and substituting the result into 



(4.3), we find 



M i2 (A) 



4CiIm(Ai 



(A - Ai)(|Ci| 2 e 2 ^i + 4e 2ie iIm(Ai) 2 ) 



Equation (3.27) yields 



16CiIm(Ai) 2 



|Ci| 2 e 2i °i +4e 2 ^Im(Ai) 2 



where Q(y, t) is defined by (3.29). Letting Ai = ia and C\ = 2ae x ° , where a > and 
io £ R are parameters, we find a real- valued solution given by 



(4.5) 



Q(y,t) 



4 arctan ( e 



The solution u(x,t) is given implicitly in terms of Q(y,t) by (3.23). We recognize 



Q(y,t) as a one-soliton solution of the sine-Gordon equation Q ty = sin(Q). We will 
see in the next section that this is no coincidence. 
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5 Relation to sine-Gordon 



The sine-Gordon equation 
(5.1) 

admits the Lax pair 



Q Ty = sin(Q) 



(5.2) 



by + i\az4> 



2 QyY) 

^((cos(Q) - l)a 3 - sin(Q)<r 2 ) 



where 4>(y, r, A) is a 2 x 2-matrix valued eigenfunction. On the other hand, the change 



of variables (x,t) i-)- (y, r) with y = p(x,t) and r = t transforms the x-part of (3.6) 
into 

i \ i l ^^XX \ 

<p y + i\oj,(p = «i H CT10. 

zvm V m / 



This x-part coincides with that of (5.2) if we identify Q y with 



1 



u x H 

m v m 



1"< 



This is the identification made in ( |3.28 | and suggests that if u(x,t) satisfies (1.1), 
then the corresponding function Q(y, r) evolves according to an equation in the sine- 
Gordon hierarchy. In fact, the following proposition shows that Q(y, r) evolves exactly 
according to the sG equation. 



Proposition 5.1 Suppose u(x,t) satisfies the generalized sG equation (1.1). Then 
the function Q(y, r) defined by 



(5.3a) 
(5.3b) 



Q(y, t) = —u(x, t) — arcsin 



u x (x,t) 



y = x — t + 



y/\+ul{x,t)) ' 

(yl + v%.{x' , t) — l)dx' , t = t, 



satisfies the sG equation (5.1). 



Proof. In view of the conservation law (3.8), we have the relations 
d d 



(5.4) 



- = - - cos(^! 



d_ 

dx 



d_ 
dy' 



Using these relations we find that the generalized sG equation (1.1) in terms of the 
variables (y, r) takes the form 



(5.5) 



u 



yr 



1 — Uy sin(u). 
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Assume that u(x,t) satisfies (1.1). Applying d T d y to (5.3a) and using (5.5) and 
its y-derivative to replace u yyT and u yT in the resulting equation, we find 



Qyr 



-u y cos(u) — yl — sin(u). 



Replacing u with —Q — arcsin(-Uy) in this equation and simplifying, we find that the 
right-hand side equals sin(Q). □ 



6 Conservation laws 



An infinite sequence of conservation laws can be derived for a bi-Hamiltonian equation 
by means of the Lenard scheme. This approach however involves the inversion of 



the bi-Hamiltonian operators. In the case of equation (1.1) the operators 9\ and 



02 are nonlocal and difficult to invert. Instead of following this procedure we can 



derive conservation laws of (1.1) by utilizing the connection with the sG equation 



established in the previous section. Indeed, suppose that the sG equation (5.1 ) admits 
a conservation law of the form 



(6.1) 



A T = B 



Vi 



where A and B are functionals of Q(y, r) and its y-derivatives. The change of variables 



(5.3b) transforms (6.1) into 



3.2) 



(\fmA) t = [B — cos(u)y/mA) 



Writing this equation in terms of u(x, t) by means of relation (5.3a), we find a conser- 



vation law for equation (1.1). Since the conservation laws of sG can be constructed 



explicitly, this approach yields an infinite number of conservatin laws for (1.1). We 
next illustrate this method by deriving a few conservation laws explicitly. 



As noted in section [2J equation (5.1) admits the recursion operator J2J1 1 where 
the Hamiltonian operators are 

Ji = d y , J 2 = dy + dyqd^qdy, 



and q = Q v . The sG equation is given by 



5H n 



where the Hamiltonian Hq is defined by 

H = J (cos^M - 1) dy. 
Further conservation laws can be determined according to the equation 



6q 
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Since J^ 1 is given by simple integration, the conserved quantities H n for n > can be 
constructed recursively. We let H n = J A n dy and write the associated conservation 
law in the form 

Hfi ~. (^4 h )t = (-^ti)j/- 

We find the following explicit expressions for the first few conservation laws expressed 
in terms of Q: 

| (a; 1 sin Q) 2 



H :{cosQ-l) T = 
Hi : Qq}) =(-cosQ) y 



H 2 ■■ [ ^Qy ~ Qly 



-Q 2 y cosQ) y , 



!l '> ' [ a Qy <fiyQyy + Qyyy 



-Q y cos(Q) - 2Q 2 y Q yy sin(Q) + Q 2 yy cos(Q) 



We use (5.3a), (5.4), and (6.2) to rewrite these conservation laws in terms of u(x,t). 



The result for Hq and H\ is 

Hq : (cos(n) + (cosu) x — Vrn) t 



H, : 



- (d x 1 (sin u + (sin u) x ) 2 + (cos(it) + (cos u) x — Vm) cos u 



(mit x + Ma) 



2m 5 /2 



(—m (m 2 + m + m x ) — u xx ) cos(n(x, t)) + 2m 2 u x sm(u(x, t)) 

2m 5 /2 



For and the expressions are quite long, thus we only present the conserved 
quantities: 

. 4 (3m x u xx - 2m(u xx + u xxx )) 21 



H 2 
H 3 



1 



8m 3 / 2 
1 

16m 5 / 2 



+ Ma; 



dx, 



m / mr 

5(mu x + u xx ) 2 (3m xU xx 2m(u xx + u xxx ) ) ( u xx 
m 6 



+ 



(m 2 (u xxx + u xxxx ) + m (-2m x u xx - 5m x u xxx + 15ul x ) - 18ul x ) 



+ ( H u x 

m 

3 A 2 1 



nr' 



dx. 



It can be verified directly that these indeed are conservation laws for equation (1.1) 



7 Traveling-wave solutions 



For a traveling wave u(x,t) = (p(x — ct) traveling with speed c, equation (1.1) takes 
the form 



(7.1) 



(cos((p) - c)ip xx = sin(yj) + ip 2 x sin(ip). 
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We rewrite (7.1 ) as 



<Pa 



sin(<^) 



cos(9?) — C 

Multiplying this equation by tp x and integrating with respect to x, we find 



(7.2) 



1 



log(l + 92 



log |c — cos(93)| + log(6), 



where b > is an integration constant. Exponentiating and then taking the square 



of both sides of equation (7.2), we find after rearrangement 

b 2 - (c - cos(<^)) 2 
(c — cos((/?)) 2 



(7.3) 



F(<p), F(tp) 



The qualitative structure of the traveling-wave solutions of (1.1) can be obtained by 



analyzing this ordinary differential equation. It is in fact possible to solve (7.3) ex- 



plicitly, which yields an expression for ip~ l in terms of elliptic functions. We proceed 
instead by briefly describing the qualitative structure of a few different types of trav- 
eling waves. We refer to [2] for a more detailed analysis in the similar case of the 
Camassa-Holm equation. 



The qualitative structure of solutions of (7.3) is determined by the zeros and the 



poles of F((p). The sets of zeros and poles of F(ip) are given by 

{<p\ cos((/j) = c ± b} and {ip\ cos(ip) = c}, 

respectively. Since b > 0, these sets are disjoint. 

The traveling waves are parametrized by b > and c £ M. Since we are interested 
in waves such that lim x ^± 00 (p(x) £ 2-7rZ we assume that b = \c — 1\. In this case 
F((p) has double zeros at ip = 2mr, n £ Z, and these zeros give rise to solutions 
(p(x) of (7.3) which approach elements in the set 27rZ exponentially fast as x — > ±00. 



In addition to these solutions there exists a range of periodic solutions obtained for 
values of b such that b ^ \ c — 1 1 . 

Note that if u(x,t) solves equation (1.1), then so do the functions (x,t) 1— > 
±u(x, t) + 27rn, n £ Z. This corresponds to the fact that F(ip) is an even function of 
period 2ir. 

7.1 Smooth waves 

Assume that b = \c — 1| and c < — 1. In this case F(ip) has double zeros at ip = 2mr, 
n £ Z, and F(ip) > away from these zeros, see Figure [TJ Moreover, the set of 
poles of F(ip) is empty. An analysis of ( |7.3[ ) shows that corresponding to the interval 
between two consecutive zeros 2nir and 2(n + l)ir there exist two smooth solutions 
(fi(x) and f2(x) such that 

lim <fi(x) = 2mr, lim <pi(x) = 2{n + 1)tt. 

x— >— OO X— >co 

and 

lim <fi2(x) = 2{n + l)ir, lim (f2(x) = 2wk. 

x— >— 00 x— >oo 

This leads to the following result. 
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Figure 1 One period of the function F((p) for b = |c — 1| and c = —2. 

Proposition 7.1 For any c < —I and ngZ, there exists: 

• A smooth traveling-wave solution u(x, t) = ip(x — ct) of such that 
lim (f(x) = 2nir, lim (p(x) = 2(n + l)ir. 



A smooth traveling-wave solution u(x,t) = ip(x — ct) of (1.1) such that 



lim (f(x) = 2(n + l)ir, lim (p(x) = 2mr. 

x — > — oo x— >oo 

The profiles of these waves for n = and c = —2 are shown in (a) and (b) of Figure 
II 

7.2 Cusped waves 



In addition to the above smooth traveling waves, equation (1.1) also admits cusped 
traveling waves (cuspons). These arise for b = \c— 1| and — l<c<l(ifc>l, then 
Ffy) ^ for all pel, and there do not exist nonconstant solutions of ( |7.3[ )). 

Assume first that — 1 < c < 0. In this case F((p) has double zeros at ip = 2mr, 
n £ 7; and F(ip) > away from these zeros. See (b) of Figure [2] for the graph 
of the numerator of F{tp) when c = —1/10. Moreover, F(tp) has double poles at 
if = ± arccos(c) + 2nir, n G 7L. Near one of the poles <po = ± arccos(c) + 2mr, the 
behavior of F(ip) is given by 



(7.4) 



F(v) 



1 



(1 + c)(tp - ifo) 2 



+ 



1 



(p tp . 



Let ^{x) be a solution of (7.3) such that p{xq) belongs to the interval (0, arccos(c)) 
and ^Px(xo) > at some point xo- Then, as x increases above xo, the value of (p(x) 
approaches (fo = arccos(c). At some point x c > xq, we have (f(x c ) = v?o and near 



this point integration of (7.3) using (7.4) shows that the behavior of <p(x) is given by 

T - 



(7.5) \<p(x) - </> | = V2 



1 + c 



1/4 



y/\ x ~ x c\ + 0(x - x c ), 



Xr 
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Figure 2 The numerator b 2 — (c — cos^)) 2 of the function F(tp) for b = \c — 1| and 
c = — 1 (left), c = —1/10 (middle), c = 1/10 (right). The black dots on the 92-axes 
denote the values of (p at which F{ip) has poles. 



In the case when ip(x) — (fo has the same sign locally for x < x c and x > x c , we say 
that t/j(x) has a cttsp at x = x c , i.e. v 3 ^) is continuous, smooth locally on both sides 
of x c , and 

lim ip x {x) = — lim (p x (x) = ±00. 

Cusps such that p(x) has a minimum at x c are also called anticusps. As x — > 
±00, y(x) decays exponentially to 0. A similar analysis applies to the interval 
(— arccos(c), 0) and the translations of these intervals by multiples of 2ir. This proves 
the following result in the case of — 1 < c < 0: 

Proposition 7.2 For any —1 <c< 1, c ^ 0, and neZ, there exists: 

• A cusped traveling-wave solution u(x, t) = ip(x — ct) of \l.ty such that 

lim (f(x) = 2mr, m&xip(x) = 2rm + arccos(c), inf (p(x) = 2mr. 

The maximum is attained at some point x c , and near x c the behavior of (p(x) 
is given by (7.5) with tpo = 2mr + arccos(c), see (c) of Figure^ 



• An anticusped traveling-wave solution u(x,t) = ip(x — ct) of (1.1) such that 
lim <p(x) = 2mr, vaimpfx) = 2nir — arccos(c), sup ip(x) = 2nir. 

The minimum is attained at some point x c , and near x c the behavior of (p(x) is 
given by (7.5) with tp = 2nit — arccos(c) ; see (d) of Figure^ 

For c = — 1 and any n G Z, there exists: 



A cusped traveling-wave solution u(x,t) = cp(x — ct) of (1.1) such that 

lim (p(x) = 2mr, maxtptx) = (2n + l)ir, inf <p(x) = 2nir. 

The maximum is attained at some point x c , and near x c the behavior of <p(x) 
is given by 



(7.6) \ip(x) - ip \ = 12 1 / 3 |x - x c \ 1/3 + 0{\x - x c 
with ipo = (2n + 1)tt, see (e) of Figure^ 



2/3 \ 



X 
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Figure 3 Profiles of different types of traveling-wave solutions u(x, t) = ip[x — ct) of 
equation (1.1): (a), (b) Smooth traveling waves with c = —2; (c) Cusped traveling 
wave with c = —1/10; (d) Anticusped traveling wave with c = —1/10; (e) Cusped 
traveling wave with c = —1; (f) Anticusped traveling wave with c = —1; (g), (h) 
Traveling waves with c = — 1 which combine the cusped and anticusped solutions. 
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An anticusped traveling-wave solution u(x,t) = ip(x — ct) of (1.1) such that 



lim if{x) = 2mr, min <p(x) = (2n — 1)tt, sup</j(x) = 2nir. 

The minimum is attained at some point x c , and near x c the behavior of <p(x) is 
given by (7.6) with ipo = (2n — l)ir, see (f) of Figure^ 



Proof. It remains to prove the cases c = —1 and < c < 1. For < c < 1, F((p) has 
double zeros at <f> = 2nir, n £ Z, simple zeros at ip = ± arccos(2c — 1) + 2nir, n £ Z, 
and double poles at ip = ± arccos(c) + 2mr, n G Z, see Figure [2] Since 2c — 1 < c for 
< c < 1, we have arccos(c) < arccos(2c — 1). Thus, F{ip) behaves qualitatively in 
the interval [0, 2n] as follows: 

1. F(ip) > in the interval (0, arccos(c)). F((p) has a double zero at ip = and a 
double pole at ip = arccos(c). 

2. F(ip) > in the interval (arccos(c), arccos(2c — 1)). F(ip) has a simple zero at 
(p = arccos(2c — 1). 

3. F(ip) < in the interval (arccos(2c — 1), 2ir — arccos(2c — 1)). 

4. F(ip) > in the interval (2n — arccos(2c— 1), 2n — arccos(c)). F(ip) has a simple 
zero at <p> = 2tt — arccos(2c — 1) and a double pole at ip = 2tt — arccos(c). 

5. F((p) > in the interval (2tt — arccos(c), 2ir). F((p) has a double zero at (p = 2ir. 

The proof when < c < 1 follows from the above properties of the function F(ip) 
and from its 2tt periodicity. Indeed, the solutions listed in the proposition arise in the 
intervals (0, arccos(c)) and (2tt — arccos(c), 2ir), as well as in the 27r-translations of 
these intervals. The intervals (arccos(c), arccos(2c— 1)) and (2ir — arccos(2c— 1), 2tt — 
arccos(c)) give rise to periodic cusped traveling waves, which are not included in the 
proposition. 

For c = —1, F(ip) has double zeros at <p = 2nn, n £ Z, F(ip) > away from these 
zeros, and F(ip) has fourth-order poles at (2n + l)-7r, n G Z. The behavior of F{ip) 
near one of the poles (po = (2n + l)ir is given by 



It follows that the behavior of a solution (f(x) of (7.3) near a point x c at which 



ip(x c ) = is given by (7.6). The rest of the proof when c = —1 follows as in the 



previous cases. □ 



Since both the cusped and anticusped waves of Proposition |7.2| when c = — 1 
satisfy cp £ (2Z + l)ir at the cusps, it is possible to combine these waves into traveling 
waves such as those shown in (g) and (h) of Figure [3} These waves are reminiscent 
of the smooth traveling-wave solutions but have infinite slope at the points where 
(p £ ir + 27TZ. 
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We emphasize that we have only shown formally that the cusped traveling waves 
are solutions of (1.1). In particular, we have not provided an appropriate weak 



formulation of ( 1.1 ) of which these waves are solutions. However, in analogy with the 



CH equation and the form of Eq. (7.1), it is expected that such a weak formulation 
does exist. Indeed, let (p(x) be a traveling wave with a cusp located at x = x c . Then 
ip xx will have a singularity at x c . However, since <p xx in equation (7.1 ) is multiplied by 
the factor cos (<£>)— c which vanishes at x c , it is possible that the product ip xx (cos(p)— c) 
is sufficiently regular to allow for a weak formulation. Actually, a direct analog of 
this mechanism is responsible for the cancellation of the singularities for the peaked 
and cusped solutions of the Camassa-Holm equation. If this analogy with CH turns 
out to be correct, it should also be possible to construct more exotic traveling- wave 
solutions of equation ( 1.1 ), such as fractal-like waves and waves with plateaus cf. [0 



We have not investigated the stability of the cusped traveling waves of (1.1). It 



is in fact unknown even in the case of the Camassa-Holm and Degasperis-Procesi 
equations whether the cusped solutions are stable (although the peaked solutions of 
these equations are known to be stable O [131 EH] ) • 
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